Parametrised modified gravity and the CMB Bispectrum 
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We forecast the constraints on modified theories of gravity from the cosmic microwave background 
(CMB) anisotropics bispectrum that arises from correlations between lensing and the Integrated 
Sachs- Wolfe effect. In models of modified gravity the evolution of the metric potentials is generally 
altered and the contribution to the CMB bispectrum signal can differ significantly from the one 
expected in the standard cosmological model. We adopt a parametrised approach and focus on 
three different classes of models: Linder's growth index, Chameleon-type models and f{R) theories. 
We show that the constraints on the parameters of the models will significantly improve with future 
CMB bispectrum measurements. 

PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 



I. INTRODUCTION 

Cosmic acceleration is one of the major challenges 
faced by modern cosmology and understanding the very 
nature of what is sourcing it is the main focus of up- 
coming and future cosmological experiments. Several 
approaches to the phenomenon of cosmic acceleration 
have been proposed in the literature, including modi- 
fications of the laws of gravity on large scales in or- 
der to allow for self-accelerating solutions in matter- 
only universes. Well-known examples of modified the- 
ories of gravity are f{R) models or the more gen- 
eral scalar-tensor theories_[6-|J| , the Dvali-Gabadadze- 
Porrati (DGP) model and its further extensions 
such as Degravitation [12!|. In the past years several au- 
thors have analyzed constraints on modified gravity, or 
more generally departures from the cosmological stan- 
dard model, both using current datasets as well as doing 
forecasts for future surveys [T3l - [28| 

In this paper we focus on certain classes of modified 
gravity and adopt a parametrised approach to forecast 
the constraints achievable from measurements of the cos- 
mic microwave background (CMB) bispectrum from fu- 
ture experiments. 

Future high-resolution CMB maps will have the ability 
of detecting higher-order correlations in the temperature 
distribution at high significance (see e.g. and refer- 
ences therein). While the CMB anisotropy distribution 
is generally expected to be Gaussian to high accuracy, 
small non-Gaussianities could be produced in the early 
universe, during inflation (commonly referred to as pri- 
mordial non-Gaussianities, see e.g. [10]) as well as be 
sourced, at a much later epoch, by the interaction of 
CMB photons with the local universe. For instance, the 
lensing of CMB photons by dark matter structure pro- 
duces a clear non-Gaussian signal in the CMB trispec- 
trum (the Fourier transform of the four-point correlation 
function), which can be used to constrain the amplitude 
of the lensing potential. Such signal, already discovered 
by the recent ACT [3l| and SPT js^] experiments helps 
in further constraining cosmological models. 



In this paper we study the implications of another non- 
Gaussian signal expected in the CMB, i.e. the one aris- 
ing from cross-correlations between lensing and the Inte- 
grated Sachs- Wolfe effect (ISW), which affects the CMB 
bispectrum, i.e. the three-point correlation function. 

The signature of the lensing-ISW (L-ISW hereafter) 
correlations in the CMB bispectrum has already been 
discussed by several authors (see e.g. (33l - l43j ) 

Assuming the standard cosmological scenario, ACDM, 
the L-ISW bispectrum should be detected at the signif- 
icance level of 4 — 5 standard deviations. In [s^ l37l. Issj 
the possibility of constraining cosmological parameters 
through this detection has been considered; in particu- 
lar it has been shown that an accurate measurement of 
the L-ISW will help in constraining the equation of state 
and the fractional density of dark energy. Here we shall 
analyze the potential of the L-ISW signal to constrain 
modified theories of gravity. In the latter, the evolution 
of the metric potentials can generally differ significantly 
from the ACDM prediction, therefore it is natural to ex- 
pect that the L-ISW bispectrum would provide valuable 
constraints on these theories. 

The paper is organized as follows. In Section |ll] we 
present the set of modified gravity models considered for 
our analysis while in Sec. IIIII we review the L-ISW bis- 
pectrum. In Section HVl we describe the analysis method 
and in Sec. |V] we present our results. We conclude in 
Section rVll 



II. PARAMETRISED MODIFIED GRAVITY 

Many models of modified gravity have been proposed 
as alternatives to ACDM, and analysing them one by 
one is impractical. The idea behind parametrised ver- 
sions of modified gravity is exactly that of encompassing 
several models into a single framework. The parametri- 
sations that we consider for our analysis cover a fairly 
large sample of theories and allow us to draw quite gen- 
eral conclusions about the constraining power of the data 
considered. 

In our analysis we fix the background to that of the 
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ACDM model of cosmology. The latter is currently in 
very good agreement with all observables constraining 
the expansion history, and many models of modified grav- 
ity can mimic it while introducing significant modifica- 
tions at the level of perturbations. Therefore, fixing the 
background to ACDM, allows us to isolate the effects of 
departures at the level of growth of structure, where we 
expect the most significant deviations. 



A. Linder model 

In [4^ Linder introduced a simple parametrisation of 
the growth of density perturbations in the linear regime, 
via a single parameter, the growth index 7 (which we will 
denote with 7^), defined through 



g{a) = efodlnaln,^ia)-'i^-l] 



(1) 



The idea is that of capturing independently the informa- 
tion from the expansion and the growth history, respec- 
tively in rim and 7^. Since in our analysis we fix the 
background to ACDM, Jim (a) is determined by that and 
the only parameter of interest will be 7^. 

In the cosmological concordance model, ACDM, as well 
as in vanilla-type dark energy models, 7l is to good ap- 
proximation constant and equal to 7l ~ 6/11. While it 
can generally be a function of time and scale, in several 
models of modified gravity it can still be approximated 
by a number, which often differs significantly from the 
ACDM value. For instance in the braneworld gravity 
of the DGP model, 7^ « 0.68 to good approximation 
over the whole history (44| . Things are more compli- 
cated for scalar-tensor models where often the time- and 
scale-dependence of 7l cannot be neglected. However, 
7l remains a powerful trigger parameter, since any de- 
viation of it from « 6/11 would indicate a breakdown of 
the cosmological concordance model. 

As a starting point for our analysis, we assume 7l ~ 
const, and forecast constraints on this simple one param- 
eter model. 



to cases in which the coupling is a linear function of the 
scalar field, i.e. oc (3i4>. A well known example of 

the latter are f{R) theories as we discuss in Sec. Ill CI 

The free parameters of these theories are the mass scale 
of the scalar field and the couplings l3i . Since we consider 
constraints from late time cosmology, we are interested 
only in the coupling to dark matter, and therefore drop 
the index i. 

While the modifications enter through the coupling 
of the scalar field to matter, and therefore change the 
energy-momentum conservation equations, it is possible 
to keep the latter unchanged and effectively absorb the 
modifications of the evolution of perturbations in the 
Poisson and anisotropy equation. The latter are com- 
monly parametrised with two functions /i and 7, as fol- 
lows 



fi{a, k)pA 



(3) 
(4) 



where pA = pS + 3^(/9 -I- P)v is the comoving density 
perturbation. Furthermore, for Chameleon-type theories 
p and 7 are well represented by the parametrisation in- 
troduced in [45} 



p{a,k) = 



1 + /3iXl k^a' 
1 + Xj Ba^ 

1 + (32X1 k'^a' 
1 + A2 fc2as 



(5) 



where the parameters /3i can be thought of as dimen- 
sionless couplings, Xi as dimensional length-scales and s 
is determined by the time evolution of the characteristic 
lengthscale of the theory, i.e. the mass of the scalar d.o.f. 
As shown in [17], in the case of Chameleon- type theories 
the parameters {/3i, Af } are related in the following way 



/3i 



2- A 



A? 



(6) 



and 1 < s < 4, so that effectively the degrees of freedom 
are one coupling and a time-evolving lengthscale. 



B. Chameleon- type models 

Chameleon-type theories correspond to gravity plus a 
scalar degree of freedom which is conformally coupled 
to matter fields, and has therefore a profile and a mass 
which depends on the local density of matter. The com- 
mon action for such theories is 



S = I d^x^ 



+5. ( 



(2) 



where (p is the scalar d.o.f., Xi is the z*^ matter field and 
Q;i(0) is the coupling of Xi to (f>. We will limit ourselves 



C. f{R) theories 



As it becomes clear in the Einstein frame, f{H) theo- 
ries are a subclass of the models described by action ([2|) , 
corresponding to a universal fixed coupling Ui — -y/2/3 0. 
Therefore they can also be described by the parametrisa- 
tion in ([5]). It can be easily seen that the fixed coupling 
a, = ^/2/3<j) gives /3i = 4/3 and (32 = 1/2. Viable f{R) 
models that closely mimic ACDM have s ~ 4 [13], there- 
fore the number of free parameters in Eqs. ([5]) can be 
effectively reduced to one lengthscale (using ^ e.g. the 
lengthscale Ai. 
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The latter is directly related to the mass scale of the 
scalar degree of freedom introduced by these theories 
and represented by fa = df /dR, known as the scalaron. 
Specifically, Ai sets the inverse mass scale of the scalaron 
today, i.e. Ai = l/mPj^. The results in the literature are 
usually presented in terms of a parameter Bq which is 
related to Ai as follows |46| 



2iJ2A? 



(7) 



Studying this particular subclass is interesting because 
some models belonging to this category have been shown 
to be cosmologically viable and pass local tests of grav- 
ity [13. 



III. THE LENSING-ISW BISPECTRUM 

As it is well known, if the expansion of the universe 
is not matter dominated {i.e. flm 7^ 1), the time vari- 
ation of the gravitational potential provides an addi- 
tional source of CMB anisotropics; restricting to the lin- 
ear regime, this effect is known as the Integrated Sachs 
Wolfe (ISW) effect 48], given by: 



ST 



(n)| 



ISW 



J dx + (n,x). 



(8) 



where r is conformal time. This adds a secondary 
anisotropy to the primordial signal and a significant con- 
tribution from it is expected at late times, when the Uni- 
verse starts accelerating. 

Furthermore, the weak gravitational lensing by mat- 
ter density fluctuations between us and the last scatter- 
ing surface shifts the observed direction of photons. In 
practice, the temperature anisotropy measured by an ob- 
server in the direction n is actually the anisotropy in the 
direction {n + V(f>{n)), i.e. 



Sf{n) = ST{n + V(j){n)) 



(9) 



where 5T{n) is the lensed anisotropy and ST{n) is the 
unlensed one (primordial plus ISW). The deflection angle 
is written in terms of the lensing potential (j), which is a 
weighted integral of the Weyl potential $ -I- ^ over the 
line of sight: 



dxgix) [(1' + 3')(n,x)] 



(10) 



Here x is the comoving distance from the observer and 
9ix) = (x* ~ X)/{XX*) with X* the distance at the last 
scattering surface. 

The Weyl potential enters both in the ISW and lensing 
kernel, therefore the two effects are correlated and they 
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FIG. 1. Dependence of the cross temperature-lensing Cf"^ 
angular spectra (top panel) and of the reduced bispectrum 
(bottom panel) on modified gravity parameters for the differ- 
ent models considered in the analysis. The solid curves corre- 
spond to ACDM, the dotted and dashed curves to the Linder 
parametrisation with 7l = 0.645, the long-dashed curves to 
an f(R) model with Bq = 0.42, the short-dashed curves to 
a Chameleon model (Cham/3i) with fii = 1.3, Bq = 0.50 
s = 2.0, while the dotted line to a Chameleon model (Cham 
s) with ^1 = 1.3, Bo = 0.50, s = 3.3. 



contribute a non-zero third order statistics in the CMB, 
i.e. the L-ISW bispectrum. 

As usual, it is convenient to consider an expansion in 
spherical harmonics of the temperature field: 



ST 
~T 



im^im ( 



(11) 



e=2 m= 



well 



as of the lensing potential, 0(n) = 
^my^m(n). Taylor expanding Eq. ^ in the lens- 



as 

ing potential, and applying the above harmonic expan- 
sions, one obtains the following relation between the 
lensed and unlensed multipole coefficients, (to first or- 
der in in the lensing multipoles): 
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where the coefficient feie2t is given by 



, _ ,'-h{h + i) + e2{e2 + i) + e{£ + i)\ „ 

/£ife£ — I ^ I i-f 



with 



(12) 



(13) 



-11(213 



(24 + l)(2£2 + l)(24 + l) 







i2 





(14) 

Then, the angle averaged CMB bispectrum generated 
by the lensing-ISW correlation is given by 



analysis from the spectrum, L-ISW bispectrum and their 
combination in order to compare the parametrised mod- 
els of Sec. HI] to a fiducial model, chosen to reproduce a 
ACDM cosmology. We fix the cosmological parameters 
according to the WMAP 7-year data best fit [1^ and 
vary only the parameters entering the parametrisations 
described in Sec. HH Spanning over the parameter space, 
we calculate the spectrum and L-ISW bispectrum using 
MGCAMB and build the likelihoods as described in the 
following. 

Each theoretical model is then compared to the fidu- 
cial model with a simple function which assumes that 
the spectrum and bispectrum can be safely described as 
Gaussian variables [sl, IsF, 
perature anisotropy spectrum we have 



For the standard Cf^ tem- 
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^TT.th ^TTJid 



(17) 
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^1 t2 t3 
mi 7712 ^^3 



{O'linii (^i2'm2 ^ism^ ) 



= hii2i3Cjj' Cl]^ + 5 perm.. 



(15) 



where Cj'^ is the temperature (primordial plus ISW) 
power spectrum and Cj'^ is the cross temperature-lensing 
angular power spectrum, Cj* = ((/)|^a^m), which de- 
pends on the Weyl potential and its first time-derivative 
(see e.g. [S^). In deriving ((T5|) we have imphcitly as- 
sumed the statistical isotropy of the universe and have av- 
eraged the three-point correlation function (in harmonic 
space) over the orientation of triangles by mean of ro- 
tational invariance. Numerical codes evolving perturba- 
tions typically work with the reduced bispectrum, defined 



where the uncertainty ct* is given by 



(2£+l) 



C, 



TT 



(18) 



We do not include any covariance noise matrix in 
Eq. ([T5)) . effectively assuming to be cosmic- variance 
limited up to ^ = 1000, which is a good approximation 
for future Planck-like experiments. 
For the bispectrum we have 



Xb 



E 



B 



th 



B 



(19) 



Bf 



where the sum is over all possible combinations of 
(16) ^1,^2,^3 with {ii < £2 < 4), ^1 + ^2 + £3 even and we 



In Figure [T] we plot different theoretical predic- 
tions for Cj'^ and the reduced bispectrum fo^^^a^s 
computed with the publicly available code MGCAMB 
( |http : //www ■ sf u ■ ca/ -aha25/MGCA MB . html) . As it can 
be noticed, the L-ISW bispectrum is clearly sensitive to 
modifications of gravity and in principle can be used to 
put constraints on models of modified gravity. In the 
next Section we describe the analysis method that we 
have used to forecast the latter for Planck-like experi- 
ments. 



IV. FUTURE CONSTRAINTS FROM CMB: 
METHOD 

We shall estimate the potential of upcoming L-ISW 
bispectrum measurements from CMB Planck-like exper- 
iments to constrain the modified gravity theories de- 
scribed in the previous section. We perform a likelihood 



set ^max = 1000, which roughly corresponds to the max- 
imum multipole sensibility for Planck-like experiments, 
(since at higher multipoles the contamination from fore- 
ground point sources starts to be dominant). 
The uncertainty cr^^ » is given by 



riTTriTTriTT 



(20) 



where ^^2^3 is 6 for equilateral configurations (£1 — I2 — 
€3), 2 for isoscele ones (with two multipoles equal) and 1 
for the scalene ones (when all the multipoles are differ- 
ent). There is no noise covariance matrix in the C|""^. 

Once the functions are computed, we can build 
the separate likelihoods for the spectrum and bispectrum 
data respectively: 



exp 



Xlh 



(21) 
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Neglecting the correlation between spectrum and bis- 
pectrum, we can further combine them in a total likeli- 
hood as follows: 



/:, = /:,A-cxp(^-^^^^j (22) 

In the calculation of the likelihood from the CMB an- 
gular power spectrum we do not include the lensing term 
that is clearly correlated with L-ISW bispectrum. Fur- 
thermore, when combining the two likelihoods like in 
(I22p . we are neglecting correlations between spectrum 
and bispectrum data that could arise from the large scale 
ISW term. As we will see in the next Section this is a 
good approximation since the bispectrum will constrain 
modified gravity parameters with a much stronger sig- 
nificance than spectrum data alone. When computing 
the bispectrum we do not include the non-linear Rees- 
Sciama term, since that would require a modeling of non- 
linearities in modified gravity. The exclusion of the RS 
term is expected to affect our results at most by ~ 17%, 
and therefore should not change our conclusions to a sig- 
nificant level. We plan to investigate in a future work the 
non-linear RS term in the framework of modified gravity. 
Finally, since we are modeling the (primordial plus ISW) 
spectrum as a Gaussian variable, we are effectively ne- 
glecting any inflationary non-Gaussian signal; further- 
more, we ignore contributions to the bispectrum from 
the lensing-Sunyaev-Zel'dovich correlation. Both signals 
could anyway be removed exploiting their different angu- 
lar dependence (see e.g. [sij). 

For each theoretical model of SecUll while keeping the 
cosmological parameters fixed to their WMAP 7- year val- 
ues, we vary the modified gravity parameter(s) (one at a 
time for the models that have more than one parameter) , 
and compute the spectrum and the L-ISW reduced bis- 
pectrum with MGCAMB ; we then use Eq. (fTH]) to com- 
pute the L-ISW bispectrum from the reduced one. We 
also choose a fiducial model, as discussed in the follow- 
ing, and compute the corresponding spectrum and bis- 
pectrum. Once a sufficient number of spectra is calcu- 
lated, we compute the likelihood profiles and extract the 
confidence levels on the parameter of interest. 

For each parametrisation, we choose a fiducial model 
based on a set of parameters that are, for most of the 
cases, the parameters that would reduce the cosmology 
to the ACDM one. In the case of the Linder model this 
is achieved by setting 7l = 0.555 Q. For f{R) theories. 
Bo = is the value giving ^ = 1 = 7 which are the values 
of these functions in ACDM. For chameleon theories the 
choice of the fiducial model is more complicated. Let us 
start employing the dimensionless parameter Bq ([7]) in 
place of the lengthscale Ai, so that the parameters for 
these models become {Bq, s ,/3i). As a matter of fact, we 
have three free parameters, no strong theoretical reasons 
to fix two of them and a complete degeneracy among 
the parameters when trying to reproduce ACDM, i.e. if 
we fix either i?o = or /?i = 1. We therefore proceed 





fiducial 


S 

95% c.l. 


B 

95% c.l. 


C 

95% c.l. 


7L 


0.555 


+0.044 


+0.060 
-0.056 


+0.034 
-0.042 



TABLE L Constraints on 7l of Linder model from the spec- 
trum (S) and bispectrum (B) and combined (C) analyses. 



by making a somewhat arbitrary choice on the fiducial 
model, fixing /3i = 1, Bq = 0.5, and s — 2 when studying 
the constraints on (3i and /?i = 1.3, Bq = 0.5, s = 2 when 
studying the constraints at varying s. 

In the case of Linder's model we evaluate the likeli- 
hoods in the range 0.475 < 7l < 0.635, at steps of 0.002 
for values near the fiducial one and at steps of 0.01 for 
values near the boundaries. In the case of f{R) we ex- 
plore the likelihood function in the range < Bq < 0.7, 
varying Bq at steps of 0.1. In the chameleon case we use 
a step of 0.01 for /3i and of 0.2 for s. 



V. RESULTS & CONSTRAINTS 

A. Linder model 

In Table [J and Figure [2] we report the constraints on 
7l from the spectrum, the L-ISW bispectrum and the 
combined analyses. As we can see the spectrum and bis- 
pectrum data are somewhat complementary: the CMB 
bispectrum is more powerful in constraining the pa- 
rameter in the region of values lower than those of the 
fiducial one; on the contrary the temperature anisotropy 
spectrum is more efhcient for larger values. The non- 
gaussian shape of the likelihood from the temperature 
spectrum can be easily understood by the fact that even 
in the case of small ISW signal (when jl 0) the an- 
gular spectrum is different from zero and still provides 
a reasonable fit to the data. The bispectrum is, on the 
contrary, not null only if the ISW is different from zero 
and it therefore provides a much more reliable way to 
detect it. 

As we can see, spectrum data provides solely an upper 
limit for 7^, leaving it practically unconstrained on the 
lower tail. On the contrary, bispectrum data give a ~ 
5% error on 7^. When spectrum and bispectrum data 
are combined there is a substantial improvement in the 
measurement. 



B. Chameleon models 

The constraints on Chameleon models from the spec- 
trum, L-ISW bispectrum and combined analyses are re- 
ported in Table HIl and Figure [31 
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FIG. 2. Likelihood distribution function for the growth in- 
dex 7l IU from the analysis of spectrum, bispectrum and 
combined data. 





fiducial 


S 

68% c.l. 


B 

68% c.l. 


C 

68% c.l. 


Pi 


1.00 


+0.25 
-0.17 


+0.10 
-0.13 


+0.09 
-0.10 


s 


2.00 


+0.55 
-0.17 


+0.42 
-0.28 


+0.30 
-0.15 



TABLE II. Constraints at 1 standard deviation on the 
Chameleon models parameters /3i and s coming from the 
analysis of spectrum (S), bispectrum (B) and combined (C) 
datasets. 
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FIG. 3. Likelihood distribution function for Chameleon mod- 
els from the spectrum, bispectrum and the combined analyses. 
The top panel gives the likelihood of the parameter /3i when 
Bo and s are fixed to Bo = 0.5 s = 2; the bottom panel gives 
the likelihood for s when Bo and /3i are fixed to Bo = 0.5 and 
/3i = 1.3. 



As for the Linder model, the two datasets are com- 
plementary in constraining the Chameleon parameters. 
The simple temperature power spectrum is more power- 
ful in constraining values of /3i < 0.75, i.e. the lower tail, 
while the bispectrum data provide stronger constraints 
on the higher tail, where the spectrum data leave the pa- 
rameter practically unconstrained. The same behaviour 
is seen for the likelihood distribution of the s parame- 
ter. Small values of s (s < 2) can be better constrained 
by temperature spectrum data. However large values of 
s are left unconstrained from the temperature spectrum 
and, on the contrary, are significantly constrained when 
using the bispectrum. This is related, as we already ex- 
plained when discussing constraints on 7^, to the entity 
of the ISW signal in the two tails; namely, the spectrum 
looses constraining power in the parameter range where 
the ISW is suppressed and tends to zero. 



C. f{R) theories 

In Table Hill and Figured] we report the constraints on 
/(i?) models from the spectrum, L-ISW bispectrum and 



combined analyses. In this case the constraints coming 
from the bispectrum are definitely tighter than the ones 
from the temperature spectrum. Once again, this is re- 
lated to the ISW signal which is suppressed w.r.t. ACDM 
one for all the values of Bq in the range < Bq < 3/2 
(becoming null at Bq — 3/2) [53|. Current constraints 
from ISW data from CMB-galaxy correlations are of the 
order of < 0.4 [26]. As we show, the L-ISW bis- 
pectrum can clearly improve CMB constraints on these 
theories, tightening the bounds by a factor of six. 





fiducial 


S 

68% c.l. 


B 

68% c.l. 


C 

68% c.l. 


Bo 





< 0.61 


< 0.14 


< 0.10 



TABLE III. Constraints at 1 standard deviation on f{R) the- 
ory parameter Bo coming from the analysis of spectrum (S), 
bispectrum (B) and combined (C) datasets. 
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FIG. 4. Likelihood distribution function for the parameter 
So describing f{R) theories from the analysis of spectrum, 
bispectrum and combined data. 

VI. CONCLUSIONS 

High accuracy temperature maps of the CMB 
anisotropy from ongoing and future experiments will 
provide an unique opportunity to test non-Gaussianity. 
While some of the signal could be primordial, a clear 
non-Gaussian signal is expected from the correlation of 
lensing and the Integrated Sachs- Wolfe effect. This sig- 
nal provides a new test of the cosmological scenario per 
se, and could further be used to test alternatives to the 
cosmological constant in the context of cosmic acceler- 
ation. In this paper we have considered three different 
parametrisations of modified gravity and investigated the 
improvement in constraining their parameters by includ- 
ing the signal in the bispectrum coming from lensing-ISW 
correlations. We have found that in the case when all the 



cosmological parameters are fixed, the bispectrum signal 
will be extremely useful providing a significant improve- 
ment in the constraints on modified theories of gravity. 
While the forecasted constraints have been obtained with 
the assumption of the cosmological concordance model as 
the fiducial one, we believe that our results have little de- 
pendence on this choice, since current data accepts only 
relatively small deviations from the standard picture. 

Finally, while the L-ISW bispectrum signal will be pre- 
sumably detected by the Planck satellite mission at about 
four standard deviations, the CMB lensing signal will 
be detected at much higher statistical significance and 
could also provide useful constraints on modified grav- 
ity theories (see e.g. [l^l)- However CMB lensing is not 
directly sensitive to time variations in the gravitational 
potentials, which instead enters directly in the L-ISW 
signal. Measurements of the ISW signal through correla- 
tions of CMB maps with galaxy surveys already provide 
interesting constraints on the models presented here [26| , 
however only a marginal future improvement in this mea- 
surement is expected. As we have shown, the constraints 
coming from observations of the L-ISW bispectrum, be- 
ing sensitive to both the spatial gradient and the time 
variation of the Weyl potential, will be complementary 
to these other observations, improving CMB bounds on 
modified theories of gravity. 
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